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Abstract 



A hybrid censoring scheme is a mixture of Type-I and Type-II censoring 
schemes. We study the estimation of parameters of weighted exponential dis- 
i -Q tribution based on Type-II hybrid censored data. By applying EM algorithm, 

maximum likelihood estimators are evaluated. Also using Fisher infirmation 
matrix asymptotic confidence intervals are provided. By applying Markov 
Chain Monte Carlo techniques Bayes estimators, and corresponding highest 
posterior density confidence intervals of parameters are obtained. Monte Carlo 
simulations to compare the performances of the different methods is performed 
and one data set is analyzed for illustrative purposes. 
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kj . 1 Introduction 

H 

Type-I and Type-II censoring schemes are two most popular censoring schemes 

which are used in practice. They can be briefly described as follows. Suppose n 
units are put on a life test. In Type-I censoring, the test is terminated when a 
pre-determined time, T, on test has been reached, and failures after time T are not 
observed. In Type-II censoring, the test is terminated when a pre-chosen number, 
R, out of n items has failed. It is also assumed that the failed items are not replaced. 
So, in Type-I censoring scheme, the number of failures is random and in Type-II 
censoring scheme, the experimental time is random. 
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A hybrid censoring scheme is a mixture of Type-I and Type-II censoring schemes 
and it can be described as follows. Suppose n identical units are put to test. The 
test is finished when a pre-selected number R out of n items are failed, or when 
a pre-determined time T on the test has been obtained. From now on, we call 
this Type-I hybrid censoring scheme and this scheme has been used as a reliability 
acceptance test in [23]. This censoring scheme was introduced by Epstin [12], he 
also studied the life testing data under the assumption of exponential distribution 
with mean life 9. Epstein [12] proposed two-sided confidence intervals for 8 without 
any formal proof. Fairbanks et al. [13] moderated partly the proposition of Epstein 
[12] and suggested a simple set of confidence intervals. Chen and Bhattacharya 
[3] earned the exact distribution of the conditional maximum likelihood estimator 
(MLE) of 9 and implied a one-sided confidence interval. Childs et al. [5] proposed 
some simplifications of the exact distribution. From the Bayesian point of view, 
Drapper and Guttmann [8] studied the same problem, and reached a two-sided 
credible interval of the mean lifetime based on the gamma prior. Comparison of 
the different methods using Monte Carlo simulations, can be found in Gupta and 
Kundu [T7]. For some related work, one may refer to Ebrahimi |101 111] . Jeong et 
al. [18] , Childs et al. [5], Kundu [19] . Banerjee and Kundu [T|, Kundu and Pradhan 
[20] . Dube et al. [9] and the references cited there. 

One of the disadvantages of Type-I hybrid censoring scheme is that there may 
be very few failures occurring up to the pre-fixed time T. Because of this, Childs et 
al. [5] proposed a new hybrid censoring scheme known as Type-II hybrid censoring 
scheme which can be described as follows. Put n identical items on test, and then 
stop the experiment at the random time T* = max{:E^ :n ,T} , where R, and T are 
prefixed numbers and xr :u indicates the time of Rth. failure in a sample of size n. 
Under the Type-II hybrid censoring scheme, we have one of the following three types 
of observations: 

Case I: {x 1:n < • • • < x R:n } if x R:n > T. 

Case II: {xi :n < • • • < x d -. n <T < xa+i-.n} if R < d < n and x d:n <T < x d+ i :n . 
Case III: {xi :n < • • • < x n - n < T}, 

where x± :n < • • • < x R:ri denote the observed ordered failure times of the experi- 
mental units. A schematic illustration of the hybrid censoring scheme is presented 
in Figure [1] 
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Figure 1: A schematic presentation for Type-II hybrid censored scheme. 

In this article, we consider the analysis of Type-II hybrid censored lifetime data when 
the lifetime of each experimental unit follows a two-parameter weighted exponential (WE) 
distribution. This distribution was originally proposed by Gupta and Kundu [15] . The 
two-parameter WE distribution with the shape and scale parameters a > and A > 0, 
respectively, has the probability density function (pdf) as: 

f WE ( X ;a,\) = ^-\e- X:c (l-e- aXx ); * > 0- (1-1) 

a 

We denote a two-parameter WE distribution with the pdf (jl.lj) by WE(a, A) and the cor- 
responding cumulative distribution function (cdf) by Fwe{x\ a, A). 

The aim of this article is two fold. First, we try to earn the MLE's of the unknown 
parameters. It is observed that the maximum likelihood estimators can be obtained im- 
plicitly by solving two nonlinear equations, but they cannot be obtained in closed form. 
So MLE's of parameters are derived numerically. Newton-Raphson algorithm is one of the 
standard methods to determine the MLE's of the parameters. To employ the algorithm, 
second derivatives of the log-likelihood are required for all iterations. The EM algorithm is 
a very powerful tool in handling the incomplete data problem see Dempster et al. [6] and 
McLachlan and Krishnan [23]. Then we use the EM algorithm to compute the MLE's. We 
also evaluate the observed Fisher information matrix using the missing information principle 
which have been used to obtained asymptotic confidence intervals of the unknown param- 
eters. The second aim of this article is to provide the Bayes inference for the unknown 
parameters for Type-II hybrid censored data. It is observed that Bayes estimators can not 
be obtained explicitly, we provide two approximations namely Lindlcy's approximation and 
Gibbs sampling procedure. So we use the Gibbs sampling procedure to compute the Bayes 
estimators, and the HPD confidence intervals. We compare the performances of the different 
methods by Monte Carlo simulations, and for illustrative purposes we have analyzed one 
real data set. 

The rest of the article is arranged as follows. In Section 2, we provide The MLE's of 
the unknown parameters. Fisher information matrix is evaluated in Section 3. Using Lind- 
ley's approximation and Gibbs sampling we obtain Bayes estimators and HPD confidence 
intervals for the parametes in Section 4. Simulation results are presented in section 5. We 
verify our theoretical results via analyzing data set in Section 6. 



2 Maximum likelihood estimators 

In this section, wc study MLEs of the model parameters a and A for WE (a, A) distribution 
with density function: 

f(x) = ^-i-Ae- Ax (l - e- aXx ) : a, A, x > 0. 
a 

For simplicity, we apply a re-parametrization as a and f3 = aX. By this, the WE(a, A) 
distribution can be written as: 

f(x) = ^±l^e-^ B (l - e-P") : a,/3,x > 0. (2.2) 

a 

The likelihood function in Case I is given by 



n\ 

{n~R)Y 



L («' A ) = 7T i! Wn£ 1 /(s i )(l-% fl) ))(- s ) 1 (2.3) 



for Case II, 



L(a, A) = ^L-n^/fcXl - F(T))(- d ), (2.4) 

and for case III, 

L(a,A)=nr =1 /(xi), (2.5) 

where /(a;) is presented by (|2.2|) . so 

F(x) = 1 - ie-^fa + 1 - e~ 0x ). 
a 

We present likelihood functions (1231) . ([2^4]) and (j2~5j) by: 

L(a,A) = r ^— IIJ =1 /( a r i )(l-i ? (c)) (n - r >, (2.6) 

(n — r)\ 

where 

i? for Case I 

d for Case II (2.7) 

n for Case III, 



and 

XR-.n for Case I 

T for Cases II and III. 



(2.8) 
Taking the logarithm of Equation 12.61 we obtain 

l(a, \)=r ln(a + 1) - (n + r) ln(a) + r m(/3) - £ V a* + V ln(l - e - ^*) 

a ■'-^ ^-^ 

+ (n - r)(--)c + (n-r) In (a + 1 - e _/3c ), (2.9) 



a 



then the normal equations are 



M. - r n+r , J_(\^r / _ n n _^ (n-r) 



?a ~~ a+1 a r a 2 VZ^i=l A ! ^ V"- ' W T a+l-e" C2 10) 



d0 - j3 



§ - ^(El=i x* + (n - r)c) + E[ =1 fz^r + (n - r)^-_ 



Maximum likelihood estimators can be secured by solving these equations, but they 
cannot be expressed explicitly. So we use EM algorithm to compute them. The advantage 
of this method is that it is convergence for any initial value fast enough. 



2.1 EM algorithm 

The EM algorithm, originally proposed by Dempster et al. [6], is a very powerful tool for 

handling the incomplete data problem. 

Let us symbolize the observed and the censored data by X = (Xi :n , ■ ■ ■ ,X r:n ) and Z = 

(Z\, • • • , Z n _ r ), respectively. Here for a given r, (Z±, • ■ • , Z n — r ) are not observable. The 

censored data vector Z can be thought of as missing data. The combination of W — (X, Z) 

forms the whole data set. In next we follow the method Kundu and Pradhan |20j for missing 

data introducing. 

If we denote the log-likelihood function of the uncensored data set by 



l c (a, P)=n ln(a + 1) - 2n ln(a) + n ln(/3) - - [ ^ X r + ^ Z, 



,i=i t=i 



\i=l »=1 / 

For the E-step of the EM algorithm, one needs to compute the pseudo log- likelihood function 
as l 8 (a,0) = E(l c (a,(3\X)). Therefor, 

l s {a,p) =nIn(a + l)-2n]n(o)+nln()8)-- (E Xi ) +zZ ]n ^ 1 ~ e ~^ 

\i=l / i=\ 

--(n-r)A(c;a,p) + (n - r)B(c;a,f3), 
a 

where 1 

A[c' ) a,P}=E(Zi\Z i >c) and B(c,a,P) = £[ln(l - e~ ?Zi )\Zi > c], 

and they are obtained in Appendix A. 

Now the M-step includes the maximization of the pseudo log-likelihood function 12.111 
Therefore, if at the kth stage, the estimation of (a,/3) is (&k,$k), then (a,k+i,$k+i) can be 
obtained by maximizing 

g(a,/3) =nln(a+l)-2nln(a)+nln(/3)- - \J2 Xi ) + E 11 ^ 1 ~ e_fe ) 



— (n - r)A(c; a k , Pk) + (n - r)B(c; a k ,Pk) 
a 



(2.12) 



Note that the maximization of 12.121 can be earned quite effectively by the similar method 
proposed by Gupta and Kundu [16] , First, (3 k +i can be obtain by solving a fixed-point type 
equation 

h(p) = P- 



The function h(j3) is defined 

h(/3) = n 
where 



B 

W) 



-Sr 



-/3xi 



-i -1 



-@Xi 



B = ^x l + (n- r)A(c, a k ,J3 k ) 



and 



a(/3) 



y/{/3B -2n) 2 +4nf3B + (f3B - 2n) 



2n 



One can follow iteration method. Once /3 k +i is determined, a k +i can be evaluated as 
a k+ i = a(/3fc+i). 

For the estimation of A, we can use the invariance property maximum likelihood esti- 
mators and obtain A as follow: 

A=|. 

a 

3 Fisher Information matrices 

One of the advantages of using EM algorithm is that presents a measure of information in 
censored data through the missing information principle. Louis [22] improved a procedure 
for extracting the observed information matrix. In this section, we display the observed 
Fisher information matrix by using the missing value principles of Louis |22j . The observed 
Fisher information matrix can be used to build the asymptotic confidence intervals. 
Using the notations: 9 = (a,/3), X=observed data, W=complete data, ix(0)=observed 
information, ivi/(#)=complctc information and /^|x($) = mfs s hrg information, follow the 
relation to 

lx(0) = lw(0) - l w \x(0), (3-13) 

to evaluate lx{0). 

Complete information and the missing information are given respectively as: 



and 



I w (0) = -E 

Iw\x{0) = -{n-r)E 



d 2 L c (W;6) 



d6 2 
d 2 lnfz(z\y, 

de 2 



(3.14) 



As the dimension of 9 is 2, Ix{9) an d Iw\x{6) are both of the order 2x2. 

The elements of matrix lyy (9) for complete data set are presented in Gupta and Kundu 
}15) . They re-parametrized WE(a, A) distribution as A and /3 = aX. 
We report I\y{9) which have been evaluated by them here as: 



lw{0) 



0-21 0-22 



where 
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(/3 


+ A) : 


A 2 




= a 2 i 
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n 



nA(/3 + A)A 



(/3 + A) 2 /3 2 



/3 4 



in which A 



1 {Hl-y)f{l-y)i 



■dy. 



On the other hand, with the above re-parametrization and by using f|3. 14[) . one can easily 
verify 

&n(c;a,/?) b 12 (c;a,P) 

b 2 i (c; a, /3) 6 2 2(c;a,/3) 



where 

6n(c;a,/3) 



Av|x(0) = (n-r) 

2/3 



1 



6i 2 (c;a,/3) = 6 2 i(c;a,/3) = 



6 22 (c;a,,0) 



(/3 + A) 2 A 3 (f + l-e-^) A 4 (f + 1 
1 1 



«* 



ce 



-p c \ 



(/3 + A) 2 A 2 (f + l-e-^ c ) 



A 2 (f + 1 



-/3c)2 



c 2 e-^ c 



(* + *■ 



-/3c\2 



(/3 + A)B 



((3 + A) 2 | + l - e -/3c 



(f + 1 



-I3c)2 



03 e -Ac(| + I 



-/3c)' 



which B = £_._,. (1 " (1 ^» 2(1 -^ ? 



in which i( = jj_ e ., t """'- '" y -~ " dy. 

Now, Ix{9) can be computed by (|3.13[) . The asymptotic variancc-covariancc matrix 
of 9 can be obtained by inverting Ix{9). We use this matrix to secure the asymptotic 
confidence intervals for A and /3. To obtain the asymptotic confidence interval for a, we use 
the non-parametric bootstrap method [25j . 



4 Bayes Estimators and Confidence Intervals 

In this section, we study Bayes estimators for parameters a and A under symmetric loss 
functions. A very well known symmetric loss function is the squared error which is defined 
as: L(f(fj,), f(fi)) = (/(/«) — /(/i)) 2 , with f(fj,) being an estimate of /(/i). Here /(/x) denotes 
some function of [i. Bayes estimators, say /.Bayes (/■*); is evaluated by the posterior mean of 
/Ou). 



Let x = (xi: n ,--- ,x r:n ) be an observed sample from the hybrid censoring scheme, 
drawn from a WE(a, A) distribution. We apply rc-paramctrization as a and (3 = aX. So 
the likelihood function becomes 



v R\ x ) CC (« + 1 )'' or -| (E[ =1 x i+ a{n-r)) / + 1 

and log-likelihood function: 



L(a,/3|x) a Ar__^ re -^ =1 «™-^( Q + i _ e -^)—n[ =1 (l - e"^*), 



/(a, /3|x) = r log(a + 1) - (n + r) log a + r log(/3) Y^ a;, + cln - r) 

a \ti J 

r 

+ {n-r) log(a + 1 - e"' 3c ) + ^ log(l - e~ Pxi ). (4.15) 

4=1 

It is assumed that /3 and a have the following independent gamma priors: 

7^ 1 (/3)oc/r' 2 - 1 e-' 9 "^ /? > 0, 

^(aJocQ'" 4- ^"™ 3 , a>0. 

So, the joint prior distribution of a and /? is of the form 

n(a, P) cc a^-ig-atos^-ig-^i ; a > 0, /3 > 0, wi > 0, w 2 > 0, w 3 > 0, w 4 > 0. 

Then the posterior distribution a and j3 can be written as 

7T(o, /3|x) = 1^4-n-r-l^+r-l^ + jw-mts^^^-^EL, x t +c(n-r)) 

x (a + 1 - e-' 3c )™- r n[ =1 (l - e- fe ) (4.16) 

where 



OO /-00 



k= I I a u ' i ~ n ~ r ~ 1 ff W2+r ~ 1 (a + l) r e~ aW3 e~P Wl e~^(^= lXi+c ( n ~ r ^ 



JO 

x(a + 1 - e- /3c )"- r n[ =1 (l - e- pxi )dadp. 

Now the Bayes estimators of a and /3 under the squared error loss function L are respectively 
obtained as: 

-i /"OO />00 

a B ayes = E[a\x] = - / / ^4-n-r^+r-l (q , + ^rg-a^g-^i 

fc Jo Jo 
xe -#(£r=i *i+c("-r))( a + ! _ e- /3c )"- T Tir =1 (l - e -P x ')dadp, 
and 

i /-OO /-00 

&*«v« = E[f3\x] = - / / a^""-'- 1 /3- 2 +''(a + iy e -^ e -^ 
k Jo Jo 



Since A is a function of a and /3, then one can obtain the posterior density function of A 
and so the Bayes estimator of A under the squared error loss function L as: 

XBayes = E[X\x] = - / / ^4+^2-n-l^+r-l^ + u y e -u(w 3 +\ Wl ) 

k Jo Jo 

xe - A (^=i x ' +c ("- r ))( u + l - e -A,M! )"- r n; =1 (l - e" Au:r ')^^. 

As these estimators can not be evaluated explicitly, so we adopt two different procedures to 
approximate them: 

• Lindley approximation, 

• MCMC method. 

4.1 Lindley approximation method 

In previous section, based on Type-II hybrid censored scheme we obtained the Bayes esti- 
mators of a, (3 and A against squared error loss function L. It is easily observed that theses 
estimators have not explicit closed forms. For these evaluation, numerical techniques are 
required. One of the most numerical techniques is Lindley's method (see [21]), that for these 
estimators can be describe as follows. In general, Bayes estimator of u(a, (3) as a function 
of a and j3 is identified: 



m = 



j °° / °° u(a, p)e l ( a ^ + P( a ^ dad/3 
/ °° /" e l (»-^)+P( a ^ dad/3 ' 



where l(a,(3\x) is log-likelihood function (defined bv !4.15p and p(a, (3) = \ogn(a,(3). 
By the Lindley's method /(x) can be approximated as: 

I(x) = u(a, (3) + -[{u aa + 2u a p a )a aa + (u Pa + 2u /3 p a )a l3a + (u afs + 2u a pp)a a0 

+ (W/3/3 + 2upf)p)app] + -[{u a & aa + UpCT a p)(l aaa a aa + l a f3a&aP + haa&pa 
+lpPa^/3p) + (u a bp a + Up<Jpp){lp aa Vaa + l a f3p&ap + lpap&/3a + ^/3/3/3 673/3)] , 

where a and (3 are the MLE's of a and (3 respectively. Also, u aa is the second derivative of 
the function u(a, (3) with the respect to a and u aa valued of u aa at (a, f3). Other expressions 
can be calculated with following definitions: 



aa z 



d 2 l 

dad/3 



2fiA 



a=&,f)=p (q, + 1)2 t a 2 a 3 (d + l-e-^ c ) 

c(n — r)e~P c 
{a + 1 - e~P c ) 



A c(n - r)e~P c 



^a.0a. — 


dH 

dadflda 


S 3 / 


-2A 2c(n - r)e" /§c 


T an 


=a, 


2r 2(n + r) | 6/3A | 2(n - r) 


""" da 3 ° 


'=" (a + 1)3 a 3 ' d4 ' (d + 1 - e-^) 3 ' 


where A = J^l-i x i + ( n — r ) c i 


f d 3 Z 


f 9 3 / 


c 2 (n-r)e-^ c (d + l + e-^ c ) 


" pp dad/3 2 


<*=«>/<=£ ip " p 9/39a9/3 


«=&,£=£ ~ (d + 1 - e-/3c)3 


f d 2 l 


—r c 2 (n — r)(a + l)e * c ^ xfe ^ 

a=a,/3=/J - ^ 2 (fi + 1 _ e -/3c)2 ^ (1 - e"^) 2 ' 


f dH 


c 2 (n-r)e~^ c (d + l + e-^ c ) 


lpp " d^da 


a~&,f>=!>- (a + l- e -/9c)3 


j d 3 l 

he? - w 


2r c 3 (n-r)(a+l)e-^ c (d + l + e^ c ) Asfe- fe (l + e- fe ) 

a=a,/3=/3 " ^ 3 ■+ ( A + 1 _ e -/§c)3 'iL (J _ e -£x«)3 


dp 


W-z — l 


« dp 

Pa = -£- 

oa 


Wi — 1 
a 


and we hg 


eve: 

















With the above defined expressions, we obtain the approximation Bayes estimators. 
Also we have: 

u(a, /3) = a, u a = 1, u aa =up = upp = u a/3 = up a = 0, 

the Bayes estimator of a under the squared error loss function L becomes 

OiBayes = & + -[2p a (J aa + 2pp(J a p + U aa l aaa + 3a aa <J a plp aa + 2(J a pl a pp 

-\-o aa applpp a + a a popplppp\. 
Proceeding similarly, the Bayes estimator of j3 under L is given by 

(u(a, (3) = 0, u /3 = 1, u aa = u a = upp = u a p = up a = 0), 

^Bayes = P + ^[2/5 Q (T Q/ 3 + 2ppOpp + 2a a/3 lp aa + 2 <J a p appl a pp 
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+°~PplpPP + ^aa^Ua + ^PP^aa^Paa + <J a p&Pplppa]- 

Finally the Bayes estimator of A under L is given by 

a a or cr a 1 

Afiayes = — + Tpaa + 2u a p a )(7 aa + (up a + 2upp a )(Jp a + (u a p + 2ll a pp)a a p 

a I 

+2upf)papp] + -[(U a a a a + Upa a p){laaaCTaa + lupuVaP + Ipau&Pa 
+lpp a O-pp) + (u a Crp a + Up(Jpp)(lpaaVaa + laPP&ap + lpapOp a + Ippp&pp)]- 

The approximate Bayes estimators of a, j3 and A can be obtained using Lindley ap- 
proximation, but it is not possible to construct highest posterior density (HPD) confidence 
intervals using this method. Therefore, we suggest the following Markov Chain Monte Carlo 
(MCMC) method to generate samples from the posterior density function, and in turn to 
obtain the Bayes estimators, and HPD confidence intervals. 

4.2 Gibbs sampling 

Here we study the Gibbs sampling method to draw samples from the posterior density 
function and then compute the Bayes estimators and HPD confidence intervals oi a, j3 and 
A under the squared errors loss function. 

Let x = (xi; n , • ■ ■ , x r:n ) be an observed sample from the hybrid censoring scheme, drawn 
from a WE(a, A) distribution. From (|4.16j) . we can write the joint posterior density function 
of a and /3 given x as: 

n{a,l3\x) oc a Wi - n ~ r ~ 1 t3 W2+r - l {a + iy e -<*ws e -p W i J e -§(£r = i*«+e(n-r)) 

x(a + 1 - e-^ c )"- r n[ =1 (l - e""")} . ( 4 - 17 ) 

by this, the posterior density function of /3 given a and x is 

7r(/3|a,x) oc j 8™"+»- 1 e H s ( tl ' 1+ £ £!"-i*«+<'(»-'-))( a + 1 - e _ * ! ) T, ~ rn <=i( 1 - e~^). 



Theorem 4.1 The conditional distribution of ft given a and x is log-concave. 

Proof 4.1 See Appendix, part B. 

By (|4.17p , the posterior density function of a given /3 and x is 

7r(a|/?,x) oc a w *-"- r - 1 (a + iy e - aW3 e -«&<=i x * +c{n - r) )(a + 1 - e^T - *"- (4-18) 
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Theorem 4.2 The conditional distribution of a given /3 and x has a finite maximum point. 



Proof 4.2 See Appendix, part C. 

Corollary 4.1 With the help of the acceptance rejection principle (see Devroye [W for de- 
tails) and the previous theorem, the generation from {4-18ty can be performed using the WE 
generator. 

Now we use theorems 14.21 and 14.21 and pursue the idea of Geman and Geman [Tl] , and 
suggest the following scheme. 

• Step 1) Take some initial value of a and /3, such as ag and (3$. 

• Step 2) Generate ct!j+i and ft+i from 7r(a|/3j,x) and 7r(/3|aj,x). 

• Step 3) Repeat Step 2, N times. 

• Step 4) Obtain Baycs estimators of a and (3 with respect to a squared error loss 
function: 

1 N . 1 N 

aBayes = _ £ «* and [3 Bayes = ]T P* 

1 »=Mi+l i=M 2 +l 

where Mi and Mi are the burn-in periods in generating of on and ft respectively. 

• Step 5) Obtain the HPD confidence interval of a: Order a±, ■ ■ ■ , aM ± as am < • • ■ < 
otiMi) an< ^ construct all the 100(1 — rf)% confidence intervals of a, as: 

("(1), a([Mi(l-»j)]))j ' • ' - («([M 1J? ])> "(Mi)), 

where [M] symbolizes the largest integer less than or equal to M. The HPD confidence 
interval of a is the shortest length interval. Similarly, we can construct a 100(1 — rj)% 
HPD confidence interval of /3. 

Finally, using the idea of Chen and Shao [4] , we can compute the estimation and HPD 
confidence interval for A. 

5 Numerical Experiments 

In this section, we carry out a simulation study to compare the performance of MLE's and 
Bayes estimators. In all the cases a = 2.5 and A = 3 arc taken. We estimate the unknown 
parameters using the MLE, Bayes estimators obtained by Lindlcy's approximations and 
also Bayes estimators obtained by using MCMC technique. We compare the performances 
of different estimators with MSE. We also obtain the average length of the asymptotic 
confidence intervals and the HPD confidence intervals. 

For computing the Bayes estimators, it is assumed that j3 and a have Gamma(i«2,wi) 
and Gamma(iC4, 103) priors, respectively. Moreover we use the non- informative priors of 
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Table 1: Average estimators, corresponding MSE and average confidence (asymptotic or 
HPD for Gibbs) length for N=40 , T=l 





R=25 


R=30 


R=35 


MLE 


2.978(4.623)3.013 
2.575(0.029)0.945 


3.014(1.836)2.993 
2.574(0.028)0.985 


2.909(0.008)2.999 
2.578(0.027)0.814 


Bayes(Lindley) 


2.863(0.019) 
2.573(0.005) 


2.968(0.019) 
2.570(0.005) 


2.863(0.019) 
2.574(0.005) 


Bayes(Gibbs) 


2.987(0.246)1.829 
2.472(0.055)0.757 


2.979(0.179)1.516 
2.482(0.053)0.751 


2.944(0.159)1.451 
2.489(0.053)0.759 



Table 2: Average estimators, corresponding MSE and average confidence (asymptotic or 
HPD for Gibbs) length for N=40 , T=2 





R=25 


R=30 


R=35 


MLE 


3.192(0.037)2.865 
2.624(0.046)0.972 


3.145(0.021)2.899 
2.591(0.033)0.891 


2.908(0.008)2.889 
2.585(0.032)0.894 


Bayes(Lindley) 


2.968(0.001) 
2.621(0.015) 


2.968(0.001) 
2.576(0.006) 


2.963(0.001) 
2.581(0.006) 


Bayes(Gibbs) 


3.148(0.216)1.544 
2.479(0.054)0.0752 


3.004(0.166)1.217 
2.482(0.054)0.748 


2.956(0.146)1.354 
2.481(0.053)0.758 



both /3 and a, by considering u>i = ?«2 = W3 = wi — 0. The Bayes estimators are computed 
under the squared error loss function and with respect to the above non-informative priors. 

The simulation is performed for different choices of n, R, T values. We replicate the 
procedure for 1000 times and report the average estimators, the MSE's, the average asymp- 
totic confidence intervals length and the average HPD confidence intervals length from the 
MCMC technique. The results are reported in Table 1-4. The first and second rows are 
parameter estimators of A and a, respectively. 

From Tables 1-4, it is observed that for fixed N and T as R increases, the MSE de- 
crease. The performances of the MLE's and Bayes estimators are very similar in all aspects. 
The average HPD confidence lengths are smaller than the average asymptotic lengths in 
all the cases considered. Finally it should be mentioned that Bayes estimators are most 
computationally expensive followed by MLE's. 

6 Data Analysis 

In this section, we demonstrate one data set for illustrative purposes. It has been studied 
by Gupta and Kundu [15] that the WE(a, A) distribution can be used quite to analyze them 
and MLE's of a and A are 1.6232 and 0.0138 respectively. The data set was studied by 
Bjcrkcdal [2] and is given below: 
12 15 22 24 24 32 32 33 34 38 38 43 44 48 52 53 54 54 55 56 57 58 58 59 60 60 60 60 61 62 
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Table 3: Average estimators, corresponding MSE and average confidence (asymptotic or 
HPD for Gibbs) length for N=50 , T=l 





R=35 


R=40 


R=45 


MLE 


3.287(0.083)2.669 
2.556(0.019)0.731 


3.079(0.006)2.667 
2.555(0.018)0.763 


3.041(0.002)2.735 
2.463(0.018)0.784 


Bayes(Lindley) 


2.979(0.000) 
2.554(0.003) 


2.984(0.000) 
2.552(0.003) 


2.981(0.000) 
2.561(0.003) 


Bayes(Gibbs) 


2.958(0.269)1.661 
2.484(0.054)0.759 


3.066(0.247)1.643 
2.485(0.054)0.757 


2.966(0.234)1.318 
2.485(0.052)0.757 



Table 4: Average estimators, corresponding MSE and average confidence (asymptotic or 
HPD for Gibbs) length for N=50 , T=2 





R=35 


R=40 


R=45 


MLE 


3.282(0.079)2.587 
2.557(0.021)0.751 


2.921(0.006)2.579 
2.558(0.022)0.745 


2.930(0.005)2.577 
2.556(0.018)0.783 


Bayes(Lindley) 


2.973(0.000) 
2.556(0.003) 


2.972(0.000) 
2.555(0.003) 


2.964(0.001) 
2.553(0.002) 


Bayes(Gibbs) 


3.194(0.214)1.418 
2.483(0.053)0.754 


3.083(0.155)1.349 
2.495(0.053)0.754 


3.001(0.045)0.504 
2.493(0.054)0.759 



63 65 65 67 68 70 70 72 73 75 76 76 81 83 84 85 87 91 95 96 98 99 109 110 121 127 129 131 
143 146 146 175 175 211 233 258 258 263 297 341 341 376. 
We use them and create the following two sampling schemes: 

Scheme 1: T = 300, R = 60, 

Scheme 2: T = 250, R = 65, 

Now for scheme 1, MLE of j3, a and A are (0.0239, 1.7715, 0.0135) and Bayes estimators with 
assumed non- informative priors, i.e., w\ = W2 = Ws = Wi = with Lindley approximation 
and Gibbs sampling method are (0.0198, 1.5019, 0.0147) and (0.0256, 2.0372, 0.0138) respec- 
tively. The 95% confidence intervals based on MLE and Bayes estimators of /3, a and A 
are 

{(0, 0.0664), (0.0073, 4.3224), (0.0075, 0.0195)} 



and 



{(0.0161, 0.0350), (1.1000, 2.9996), (0.0053, 0.0256)} 



respectively. 

For scheme 2, MLE of /3, a and A are (0.0255, 1.9390, 0.0132) Bayes estimators with assumed 
non-informative priors, i.e., w± = u>2 = W3 = w± = with Lindley approximation and Gibbs 
sampling method are (0.0223, 1.7395,0.0142) and (0.0254,2.0835,0.0131) respectively. The 
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Figure 1: Estimated density functions with informative priors for scheme l(left) and scheme 2(right). 

95% confidence intervals based on MLE and Bayes estimators of /3, a and A are 

{(0, 0.0677), (0.0024, 4.5216), (0.0075, 0.0189))} 

and 

{(0.0170, 0.0340), (1.2242, 2.9353), (0.0085, 0.0222)} 

respectively. 

Because we see the effect of the hyper parameters on the Bayes estimators and also 
on confidence intervals, we take the following informative priors w\ =3, w^ = 1.5, wy, = 
0.01, W 4 = 1. 

Based on this, for scheme 1, Bayes estimators of (3, a and A with Lindley approximation and 
Gibbs sampling method are (0.0233, 1.8233, 0.0142) and (0.0255, 1.8205, 0.0150) respectively. 
The 95% confidence intervals based on Bayes estimators of j3, a and A are 

{(0.0162, 0.0350), (1.0042, 2.6231), (0.0062, 0.0269)}. 

For scheme 2, Bayes estimators of /3, a and A with Lindley approximation and Gibbs 
sampling method are (0.0231, 1.8162,0.0141) and (0.0246,2.0799,0.0127) respectively. The 
95% confidence intervals based on Bayes estimators of j3, a and A are 

{(0.0170, 0.0322), (1.2043, 2.9135), (0.0060, 0.0217)}. 

We plot all the different estimated density functions with non-informative priors and 
informative priors in Figure 1 and Figure 2. 

Comparing the two schemes with informative and non-informative priors, it is observed 
that for scheme 1, estimators have smaller standard errors than scheme 2, as expected. Also 
it is clear that the Bayes estimators depend on the hyper parameters. Because the HPD 
confidence intervals based on informative priors are slightly smaller than corresponding 
length of HPD confidence intervals based on non-informative priors, therefore the prior 
informative should be used if they are available. 
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r lgure 2: Estimated density functions with non-informative priors for scheme l(left) and scheme 2(right). 



7 Conclusion 

In this article, we have studied the classical and Bayes inference procedure for the Type-II 
hybrid censored WE(a, A) distribution. We provide the maximum likelihood estimators and 
it is observed that the maximum likelihood estimators of the unknown parameters can not 
be obtained in the closed form and we suggest the EM algorithm to compute them, we 
also earn the Bayes estimators of the unknown parameters and show that they can not be 
obtained in explicit forms, and we have proposed two approximation methods to earn them. 
We have compared the performance of the different methods by Monte Carlo simulations, 
and it is observed that the performance of quite satisfactory. 
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Appendix A 



Theorem 7.1 Given Xn\ = X(±),--- ,Xr r \ = Xr r ), the conditional distribution of Zi for 



(r) <"(r 



where 



fz\x(Zi\X(i) =£(!),••• ,X {i) =X(i)) = fz\x{Zi\X(i) = X{i)) = 1 '_ p 7^) 






^jj(a;) = 1 - -e-° x [a + 1 - e"' 3 *] 
a 



Proof 7.1 T/ie proof can be obtained similarly as in Ng et al. (2002). 
Note that using Theorem we can write 

A(c]a,P) =E[Z l \Z. l >c] = r^s^i^ie-l^l-e-^da: (if = e -« c (a + 1 - e^ c )\ 

<a+l)a roo _„, a poo — « j ( + j9(a+l) 

= ^F-Jl c Me dM ~ A-/3(a+i) J ^CH-D ^e dv (put u = §x & » = -4-i: 
= ^(-e^Cti + l)]f c " T^TT) ("^)(« + l)]%a c 






RT/3 

S a 

K/3 



g'= /) 3(a+l)c 



( a + l)(£ c + 1) - ^__(^i±ii£ + i) 

^ / v a ' a + 1 ^ a ' 



(a + l)(/3c + a) - ^rj- (/3(a + l)c + a) 



(a + l)(^ + a)- e "' C( ^ (a+1)+tt) 



a+1 



and about B(c; a, /3), we have: 

f c 
B(c;a,/3) = E[\n(l - e~^ Zi )\Zi > c] = 

a + 1 f 1 



Q + l 



- ln(l - e- px )e-« x {l - er px )dx 



aK J i 



a K 
((l-u)° _1 ln(u) 
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^ [l^ 1 - e ~ PC ) 3 (-~ + l)hypergeom 3 . 2 ([3, 3, -1 + 2], [4,4], 1 - e^ c ) 
aK [18 a a 

-ihypergeom 2 x ([2, -- + 1], [3], 1 - e~^) 
4 ' a 

+i ln(x)hypergeom 2)1 ([2, -i + 1], [3], 1 - e ^ c )) (1 - e^ c ) 2 

-4(-- + l)hypergeom 3 2 ([3, 3, -- + 2], [4,4], 1) 
18 a a 

+ -hypergeom 21 ([2, + 1], [3], 1) 

4 a 

where hypergeom(.) is Generalized hypcrgeometric function. This function is also known 
as Barnes's extended hypergeometric function. The definition of F Piq (n,d,£) is: 

p ' q[ ' ' 4 j " ^ r(* + i)nf =1 r(d, + fc)r-i(*) ' 

where n = [ni, • • • ,n p ], p is the number of operands of n, d = [di, • • ■ ,d q ] and q is the 
number of operands of d. Generalized hypcrgeometric function is quickly evaluated and 
readily available in standard software such as Maple. 

Appendix B 

The conditional density of (3 given a and x is 

This function is log-concave because we have 



d 2 log(7r(/?|a, x)) w 2 + r - 1 (n - r)c 2 e~^ c {a + 1) ^ xf e" fe 



£<r 



<0. 



a/3 2 ,3 2 (a + l-e"^) 2 ^(i- e -fe)2 

Therefore, the result follows. 

Appendix C 

The conditional distribution of a given (3 and x is 

7r(a|/3,x) oc d^-^'-'la + l) r e - awa e--^^ x ' +c{n - r)) (a + 1 - e~ Pc ) n - r . 

In this function, we have 7r(oo|/3,x) — » 7r(0|/3,x) = and 7r(a|/3,x) > Va, now it is enough 
that prove ir(a\/3, x) is bounded. With simple calculation we see that Va this function is less 
than the gamma function and the gamma function is a bounded function, so this function 
is bounded. Therefore 7r(a|/3,x) has a finite maximum point. 
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